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Abstract. 

The raise and peel model is a one-dimensional stochastic model of a fluctuating 
interface with nonlocal interactions. This is an interesting physical model, in this 
paper we review its properties. It's phase diagram has a massive phase and a gapless 
phase with varying critical exponents. At the phase transition point, the model exhibits 
conformal invariance which is a space-time symmetry. Also at this point the model has 
several other facets which are the connections to associative algebras, two-dimensional 
fully packed loop models and combinatorics. 
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1. Introduction 

The raise and peel model (RPM), reviewed in this paper, is a one-dimensional 
adsorption-desorption model of a fluctuating interface [H [2]. The interface evolves 
following nonlocal Markovian dynamics. Originally, the model " appeared" [3] during an 
investigation of some intriguing connections between the groundstate wavefunctions of 
the XXZ quantum chain and two-dimensional dense 0(n = 1) fully packed loop models 
[HE]. For some special boundary fully packed loop (FPL) models [Ell?] are related [8] 
to Alternating Sign Matrices [9J, which are of interest in combinatorics. These magic 
connections were found by Razumov and Stroganov [T0]-|12]. 

The stationary states probability distribution functions (PDF) of the RPM are 
given precisely by the Razumov- Stroganov wavefunctions. It was later understood that 
the RPM is in fact a very interesting model on its own. Moreover, objects with a nice 
mathematical structure, and this indeed the case of Razumov-Stroganov wavefunctions, 
allow, using small systems sizes, to make conjectures for the expression of physical 
quantities for any systems sizes. 

The RPM is not only special because of the nice properties of the stationary 
states PDF but also for its dynamics. For large system sizes, the flnite-size scaling 
eigenspectrum of the Hamiltonian describing the continuous time evolution of the 
RPM is expressed in terms of characters of the Virasoro algebra [I3], hence conformal 
invariance (this implies that the dynamical critical exponent has the value z = 1). To 
our knowledge, the RPM is the flrst known example of a stochastic process having this 
space-time symmetry. As we are going to see, this symmetry has consequences for the 
physical properties of the model including the stationary states. 

By changing the ratio of the adsorption and desorption rates of the model, the 
RPM can be taken away from the Razumov-Stroganov point (which corresponds to 
the case where the rates are equal). We show the phase diagram of the model thus 
obtained. If the adsorption rates exceed the desorption rates, one gets a gapless phase 
with continuously varying critical exponents. If the desorption rates are higher than the 
adsorption ones, the system is massive (flnite correlation lengths). 

The presentation of the different facets of the RPM will be kept as elementary as 
possible, for the reader who wants to know more, we give references. 

In Section 2 we give a connection between associative algebras and stochastic 
processes. This connection is important since the structure of the RPM and some of its 
extensions are related to the Temperley-Lieb (TL) algebra [TH [15] and its extensions 
[8]. Moreover one can use two representations of the TL algebra. One in terms of RSOS 
paths which is useful for the interface model and another one in terms of up and down 
spins. In the latter representation, one gets a quantum spin chain Hamiltonian which is 
integrable and therefore the spectrum can be computed exactly using the Bethe ansatz. 
In Section 3 we give the flnite-size scaling limit of the Hamiltonian eigenspectrum. It is 
at this point that conformal invariance enters the picture. 

In Section 4 we show how the Hamiltonian expressed in terms of the generators of 
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the TL algebra acts in the representation given by hnk patterns (equivalent to RSOS 
paths). We consider the example of a Hamiltonian containing 5 generators acting in a 
space of 5 link patterns configurations. We compute the PDF of the stationary state 
(ground-state wavefunction) and show how it can be related to the 26 configurations of 
a fully packed loop model (FPL) on a rectangle or to 26 vertically symmetric alternating 
sign matrices. This phenomenon which is vahd not only for the case of five generators, 
is very interesting since it relates the stationary state of a one-dimensional stochastic 
process to an equilibrium two-dimensional system. This is the Razumov-Stroganov 
conjecture. 

All these considerations were valid for the RPM at the Razumov-Stroganov point. 
In Section 5 we describe the one-parameter dependent RPM. The parameter denoted by 
u is the ratio of the adsorption to dcsorption rates. For u not equal to one, the model 
is not anymore related to the TL algebra, few analytic calculations arc possible and the 
study is based on Monte-Carlo simulations. We give the phase diagram of the model 
which is based on the results given in the next sections. 

In Section 6 we give the properties of the stationary states, stressing what is special 
for the case u — 1. Section 7 deals with the dynamics of the model. Since the desorption 
processes are not local, avalanches occur. We briefly mention their properties. 

Some results not described in this paper are briefly presented in Section 8. We also 
list some open questions. 



2. Markovian dynamics and associative algebras. 

The continuous time evolution of a system composed by the states a — 1,2, ... ,N with 
probabilities Pa{t) is given by a master equation that can be interpreted as an imaginary 
time Schrodinger equation: 

J^Pa{t) = -Y.Ha,,P,{t), (1) 

where the Hamiltonian H is a,n N x N intensity matrix: Ha,b nonpositive (a ^ b) 

and ^aHa,b = 0. —Ha^ is the rate for the transition \h) \a). The ground-state 
wavefunction of the system |0), i?|0) = 0, gives the probabilities in the stationary state: 

|0) = $^P„|a), P„=limP„(t). (2) 

— t— >oo 

a 

The normalization factor of the unnormalized probabilities Pa is (0||0) where 

(OH^(a|, m^Q. (3) 

a 

Since H is an intensity matrix the real parts of the N eigenvalues E{k) are 
nonnegative and, if complex, the eigenvalues come in conjugate pairs. The eigenvalue 
zero is not degenerate. Conversely, if the spectrum of a matrix has these properties, 
through similarity transformations the matrix can be brought to the form of intensity 
matrices (the solutions are not unique). 
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Let us now consider an associative algebra with generators A{i) {i = 1,2, ... ,m). 
Taking products of the generators one gets the words w{r). The algebra is defined by 
giving some relations between the words. If we can choose the independent words W{r) 
such that any product of them verify the relation: 

W{r)W{s) = J2p'^'^W{q), = 1' (4) 

then the Hamiltonian 

H = J2aiu)il-Wiu)), (5) 

u 

acting in the vector space defined by the basis of independent words W{s), is an intensity 
matrix if the coefficients a{u) are nonnegative. H is acting from the left on the words 
W{s) of the vector space: W{u)W{s) (W{u) is one on the terms in ([5]) and W{s) belongs 
to the vector space). The action of the independent words W{u) on the vector space 
defined by the same words gives the regular representation of the algebra. 

If the algebra A contains a left ideajjl defined by the words I a, the Hamiltonian 
([5]) acting on this ideal gives again a stochastic process. If A has several ideals 
/i(/2(- ■ ■ (/„, H has a block triangular form. The ground-state wavefunction |0 > is 
a linear combination of the states (words) which define the vector space Ji. These 
states are also called recurrent states. The words belonging to the ideals I2, ■ ■ ■ ,In but 
not belonging to Ji do not appear in the stationary state. 

We will give a simple example. Consider the Abelian algebra given by m generators 
A{i) (z = 1, 2, . . . , L) satisfying the relations: 

A{i)^ = aA{i + iy + bA{i)A{i + l), i = 1, 2, . . . , L - 1, (6) 



A{Lf = a + bA{L), (7) 



where 



[A{{),A{j)] = 0, ^,j = l,...,L, a + b = l, a,b>0. (8) 
The model is defined by the Hamiltonian: 

H = l-jj:A{z). (9) 

^ 1=1 

H acts in the 2^ dimensional vector space of independent words which is given by 
the monomials: 1, A{i), A{i)A{j) {i j), . . .. The physical interpretation of the vector 
space is obvious. Take a one-dimensional system with L sites. Each site can be empty or 
occupied by at most one particle. In the state A{ii)A{i2) ■ ■ ■ A{ik) the site i„ is occupied 
if A{in) appears in the expression of the monomial. In this sandpile model (see Dhar's 
papers [16j for much more on this subject), the toppling mechanism is encoded in the 
algebra (I6])-(IH]) and the dynamics in Eq. ([9]). In a unit of time a particle is introduced 
with probability 1/L at any site of the lattice (see Eq. IQ). Without increasing the 
time, a given site i with more than a single particle, has the probability a (probability 

I A left I A ideal of an algebra A is a subspace of A such that aX G I a for all a e A and X Ia- 
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b) of sending two (one) of its particles to the site i + 1 (see Eq. (^). Particles at the 
site L may leave the lattice (see Eq. ([7])). The process continues up to when we have at 
most a single particle in any lattice point. 

A special class of associative algebras in which the condition (jlj) is satisfied are 
semigroup algebras. In this case, the relations (jlj) take the simple form: 

U{q) = W{s). (10) 

An important example of this kind is the Temperley-Lieb (TL) algebra with the 
generators Cj (j = 1, . . . , L — 1) satisfying: 

{q+q~^)ej, ejej±iej = Cj, CjCk = CkCj for \j-k\ > 1, (11) 



e 



with q having the special value q = exp(i7r/3). The special choice in (j5j) where we 
take a{u) = 1 and for W{u) only the generators Cj gives the exact integrable quantum 
Hamiltonian 

^=EV-e,). (12) 



For other values of q, the relations (iTTj) define an associative algebra but not a semigroup. 
The TL semigroup fjTTj) can be generalized in several ways. One can add boundary 
generators [I7]-[20] or one can generalize the group algebra itself. Examples of the 
latter are the multi-colored versions of the TL semigroup [211 E2] or the rotor model 
[23] . Physical applications of these models have not yet been studied. 

Another important example of a semigroup with possible applications to stochastic 
processes is the Brauer semigroup algebra [211 [251 126] . 

In all these cases, the unnormalized probabilities Pa which appear in the ground- 
state wavefunction |0 > are positive numbers. If these numbers are integer multiples of 
the lowest one, they might have a combinatorial meaning and the challenge is to find 
and understand them. 

From now one we will study only the applications of the TL semigroup defined by 
Eq. ( ITTi) and take the Hamiltonian ( |T2i) with L an even number {L = 2n). 

One has to specify the vector space in which the Hamiltonian (JT2j) acts. The TL 
algebra has a left ideal Ji defined by the words [151 IS] 

L/2 

W{s)J,, Jo = n 62,-1, (13) 
i=i 

where W{s) is any word of the algebra. There are 
n + l\n J 

words in this ideal. For example, if L = 4, the two independent words in this left ideal 
are: e(l)e(3) and e(2)e(l)e(3). In Section 4 we discuss in detail the action of H in this 
vector space. The raise and peel model at the Razumov-Stroganov point is defined by 
the Hamiltonian (JT2j) acting on this ideal. 
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As opposed to the Abehan algebra (jS])-® or other Abehan sandpile models, it is 
not clear at all how the TL algebra which is not Abelian can be related to a toppling 
process. 

There is another useful (reducible) representation of the TL semigroup in a spin 
basis. The matrices 

1 1 
= 2 [""i "^i+i + ""i+i " 2''^''^+^ ^ '^^""^ " ""^^^^^ ' ^^^^ 
where a^^a^ and are Pauli matrices, satisfy fill I) . The Hamiltonian given by (1121) 
becomes a spin 1/2 quantum spin chain acting in a 2^ dimensional vector space. It is 
known that this reducible representation contains all the irreducible representations of 
the TL algebra, in particular the one given by the left ideal discussed above. 



3. The raise and peel model and conformal invariance 

The quantum spin chain defined by the Hamiltonian (|T2|) and f|T5|) commutes not only 
with 

s' = \j:^-. (16) 

but with two other operators and S~ (see [27]). The three operators S*^, 5*+ and 
which are related by commutation relations define the quantum algebra Uq{sl{2)) 
for q = exp(i7r/3). This is a deformation of the usual angular momentum sl{2) algebra 
with irreducible representations of dimension 2s + 1 labeled by the spin s (integer or 
half integer). The number of scalars (s = representations) for a chain with L = 2n 
sites is given by the Catalan numbers f|T^ which coincide with the number of states in 
the left ideal given by f[T^ . It can be shown that this is not a simple coincidence. The 
spectrum of the Hamiltonian defined in this left ideal is obtained considering only the 
scalar sector of the quantum spin chain. 

The quantum spin chain is an integrable system [28]. This implies that it's spectrum 
can be computed exactly using the Bethe ansatz. If we denote by -E.^ (r = 0, 1 . . . , 2^) 
the energy levels in nondecreasing order: Eq = < Ei < E2 < ■ ■ ■, the finite-size scaling 
partition function of H is defined as follows : 

Z{q) = lim Ziiq) = lim Vg^-^"/"''% (17) 

n 

where Vg = 3^/3/2 is the sound velocity [28]. One can show [29] that Z{q) has the 
expression 

Z(g)=$:(2s+l)C(g). (18) 

s 

Here s is the spin, taking the values s = 0, 1, 2, . . ., for L even, and s = |, |, |, . . ., for 
L odd, and 

00 

n=l 
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where 

A. ^ (20, 
This imphes that for large lattice sizes, the energies are (see f llSp and fll9p ) 

E=^-^{A, + k) + o{^), (21) 

where k is an integer. One can observe that the finite-size scaling spectrum of the 
Hamiltonian is expressed in terms of characters of a Virasoro algebra [13j with a central 
charge c = (the ground-state has energy zero without finite-size corrections). Hence 
the system is conformally invariant. This implies a complete knowledge of the finite-size 
scaling spectrum. In the spin zero sector, for L even, one gets for the first excitations 
(see Eqs. (EOJ-iEID): 



E^ = 3nV3/L, = vr^, . . . , (22) 



which correspond to the values = and A; = 2, 3, ... in Eq. fl2T|) . 

In a conformal theory the A^ given by Eq. fl2U]) are related to the critical exponents 
of various correlation functions [131 • Notice that 



Ai = ^, A2 = 2 (L even) 

As = 1, A^ = — (L odd). (23) 
2 2 3 

Eqs. ([22D and fl2^ will be used in the next sections. The dynamical critical 
exponent z of a stochastic process is defined by the finite-size scaling behaviour of 
the Hamiltonian spectrum: 

hm L'E = const. (24) 

L—*oo 

Comparing Eqs. fl2Tl) and (|2^ we conclude that in the RPM, z = 1. 

4. The raise and peel model: combinatorial facets 

The TL semigroup algebra ffTTl) : 

e'j = Gj, ejCjiiej = Cj, CjCk = CkCj for \j — k\ > 1, (25) 

can be understood in terms of graphs [15J. The generators Cj can be pictorially 
represented by 

ej = \"":\": - (26) 



1 2 j-1 j j+lj+2 L-1 L 



The words in the left ideal discussed in Section 2 (Eq. f|T3|) ) can be represented 
by boundary diagrams of loops or link patterns [I5]. An example of such a diagram is 
shown in Fig. 1. If the TL semigroup has L — 1 generators one takes L sites. Contour 
lines connect pairs of sites and don't intersect. 
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Figure 1. An example of a link pattern for L = 16 



The action of Cj on a link pattern of contour lines is given by placing the graph of cj 
underneath that of the link pattern one, removing the closed loops and the intermediate 
dashed line. Next, one contracts the links in composite pictures. The action of ci on 
one of the link patterns for L = 6 is given by. 



We consider the action of the Hamiltonian 

L-l 



(27) 



(28) 



=1 



on the link patterns. As an example we consider the case L = 6. In this case we have 
five words in the ideal Ii. 1) Jq = 616365, 2) 62J0) 3) c^Jq, 4) 6264J0 and 5) 636264J0. 
These words correspond to the link patterns: 



1 : ^^....^^....^^ 





(29) 
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-4/ 



(30) 



Which is indeed an intensity matrix. In the basis (!29l) the stationary state |0) of if is 
given by 

|0) = (11,5,5,4,1). (31) 

Notice that all the components of |0) are integer numbers (we have chosen the smallest 
component to equal 1). The normalization factor (0||0) (equal to the sum of the 
components of |0)) is equal to 26. Below we show that the normalization factor acquires 
an extra meaning from an enumeration problem |3] : it is equal to the partition function 
of an equilibrium statistical mechanics system in two dimensions. This observation has 
deeper consequences which we will now briefly discuss. 
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We consider a six-vertex model [20] in a (L + 1) x (L + 1) square lattice with domain 
wall boundary condition [31]. This boundary condition imposes that all the arrows on 
the vertical (horizontal) boundaries point outward (inward). In Fig. [2](a) we show the 
fixed links in this boundary condition for the case L = 6. The arrow configurations on 
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Figure 2. a) Domain wall boundary condition for the six-vertex model, b) 
Fixed arrows for the horizontally symmetric vertex configurations. The lattice size 
is (L + 1) X (L + 1) with L = 6. 



this lattice can be transformed into FPL configurations. The latter are configurations 
of paths such that every site is visited by exactly one path. This transformation is done 
by dividing the square lattice into its even and odd sublattices denoted by A and B, 
respectively. Instead of arrows, only those edges are drawn that on sublattice A point 
inward and on sublattice B point outward as in Fig. [3l We take the vertex in the upper 
left corner to belong to sublattice A. 



■> 



A 



■< 



< <- 



A 



A 



>- 



A 



A 



-L J. 



r 



B 



r 



Figure 3. FPL vertices on sublattices A and B derived from the six arrow vertices. 



The domain wall boundary conditions of the six-vertex model translates into a 
boundary condition for loops. For example, the boundary condition shown in Fig. [2]^a) 
translates into the boundary conditions given in Fig. ID^a). 

Among the general configurations of the vertex model, we now restrict ourselves to 
those which are horizontally symmetric. In the case of L = 6 the boundary condition is 
shown in Fig.[2](b). Due to the correspondence of vertices given in Fig. [3]the related FPL 
configurations will also be horizontally symmetric. In Fig. Hl^b) we show the boundary 
condition for FPL related with that of Fig. [2](b), for L = 6. 
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U 1 . 1 J 




Figure 4. a) Boundary conditions (fixed edges) for the FPL diagrams . b) Fixed 
edges for the horizontally symmetric FPL diagrams. The lattice size is (i + 1) x (L + 1) 
with L = 6. 

The above relations show that number of horizontally symmetric configurations of 
the six-vertex model in the (L + 1) x (L + 1) square lattice with domain wall boundary 
condition is the same as the number of FPL diagrams in the (L — 1) x L/2 rectangular 
lattice. The boundary condition for the FPL diagrams is special, for the case L = 6 is 
shown in Fig. [5l 



~ I ■ I ~ 

Figure 5. FPL grid corresponding to Fig. 2. 

The paths that start and end on boundary sites (we disregard the closed loops in 
the bulk of the FPL diagrams, as a consequence of the correspondence given in Fig. [3]) 
define link patterns having the same topology as the TL link patterns (l29l) . There is a 
fascinating connection between the components of the stationary state |0) (see Eq. (|3T1) ) 
and the enumeration of FPL configurations on the rectangular grid: the unnormalized 
probability of a state corresponding to a given link pattern is equal to the number of FPL 
configurations with the same link pattern l5]. This is an yet unproven conjecture. 

Take for example the stationary state fl3Tl) for L = Q. The FPL configurations on 
the 5x3 rectangle are shown in Fig. 6. Their total number is 26, which is equal to 
(0||0) and they can be categorized according to the five link patterns present in ( l30l) . 
One finds that the number of diagrams corresponding to the link pattern 1 in (l29l) is 11 
(they are printed in bold in Fig. 6), to link pattern 2 is 5, to link pattern 3 is 5, to link 
pattern 4 is 4 and to link pattern 5 is 1. 

We have illustrated for L = 6 the connection between the unnormalized 
probabilities in the stationary state of various link patterns in the RPM at the Razumov- 
Stroganov point and the number of FPL configurations with the same link pattern. 
This connection is important since it shows that one can see a "far from equilibrium" 
stationary state of an one-dimensional system as an equilibrium state of a system in two 
dimensions. 

Before we further elaborate on this point, we will shortly discuss another facet of 
the RPM: its connection with the enumeration of alternating sign matrices (ASM) (see 
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Figure 6. The 26 FPL diagrams for L = 6. The 11 diagrams corresponding to link 
pattern 1 in ([^ are printed in bold. 



[H] for an excellent introduction to this subject). 

An alternating sign matrix is a square matrix of Os, Is and —Is for which a) the 
sum of the entries in each row and each column is 1. b) the non-zero entries of each row 
and of each column alternate in sign. An example of such a matrix is 



(32) 



Alternating sign matrices can have symmetry properties. For example, the matrix 
fO 1 o\ 

1 -1 1 (33) 
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is vertically and horizontally symmetric. There is a one to one correspondence between 
the six-vertex configurations with domain-wall boundary conditions and ASM. To each 
vertex you attach the numbers 0, —1 or 1 as shown in Fig. 7. To six- vertex configurations 
with a certain symmetry correspond alternating sign matrices with the same symmetry. 



The number of six-vertex configurations on a (L + 1) x (L + 1) square with domain 
wall boundary conditions that are invariant under reflection in the horizontal symmetry 
axis is equal to that of the (L — 1) x L/2 rectangle with special boundary conditions 
like in Fig. [2](b). The total number of such configurations is known and is equal to 
the number of (L + 1) x (L + 1) horizontally symmetric ASM (equal to the number of 



Different facets of the raise and peel model 



12 



< > 



-> < ( f 







-1 



0-11 

Figure 7. Correspondence between vertices and entries of alternating sign matrices. 



vertically symmetric ASM) which is given by [32]: 

(2j + l)!(6j + 3)! 



^2n+l 



n(3j 

j=0 



2)- 



1,3,26,646, 



(4j + 2)!(4j + 3)! 
where L = 2n. The leading asymptotic terms of A^_^^ are given by 



(L-l)L 



(L - 1) In 



3^6 



lnL2 + 0(l), 



(34) 



(35) 



2 ' ' 8 144 

where Sq = In(^j^). The first term is proportional to the area of the rectangle, the 
second one is related to the boundary. This shows that, indeed, the stationary state of 
the RPM is related to a bona fide two-dimensional statistical mechanics system. 

In the next section we are going to discuss physical applications of the stochastic 
model defined by the Hamiltonian H (Eq. fl28l) ) acting on the vector space of link 
patterns (see, for L = 6, Eq. ( l29l) ). In the Appendix we show how each link pattern can 
be mapped onto a RSOS (Dick) path. Such a path can be seen as an one-dimensional 
interface separating clusters of "tiles" deposited on a substrate, from a rarefied gas of 
tiles. The action of the Hamiltonian H on the RSOS configurations gives the transition 
rates for various adsorption and desorption processes (see Eqs. f l29|) . fl30|) and Fig. IA2I 
for an example with L = 6). 



5. The raise and peel model: the phase diagram 

We consider a one-dimensional lattice with L + 1 {L = 2n) sites. An interface is formed 
by attaching at each site non- negative integer heights hi [i = 0, 1, . . . , L) which obey 
the restricted solid-on-solid (RSOS) rules: 

hi+i - hi = ±1, ho = hL = 0, hi > 0. (36) 

There are C„ = (2?T,)!/((r;, + l)(n!)^) possible configurations of the interface, that 
also correspond to the number of independent words (HM in the left ideal Ji given in ( |T3l) 
(see Sec. 4). In Fig. [T]we show a configuration for n = 8 {L = 16). Alternatively, one 
can describe the interface using slope variables Si = (/ij+i — hi^i)/2, {i = 1, L — 1). 

The dynamics of the interface is described in a transparent way in the language 
of tiles (tilted squares) which cover the area between the interface and the substrate 
{h2i = 0, h2i+i = 1, (« = 0, ...,n-l),hL = 0) (see Fig. El- 

We consider the interface separating a film of tiles deposited on the substrate from 
a rarefied gas of tiles. The evolution of the system in discrete time (Monte-Carlo steps) 
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Figure 8. A configuration of the interface with three contact points and two clusters, 
as defined in Sec.[6l for the lattice size L = 16. 




12 3 16 



Figure 9. A dcsorption event. The incoming tile at site 1 triggers, with probability 
Md, an avalanche of 5 tiles, which are shaded. All of the shaded tiles are removed in 
the desorption event. 



is given by the following rules. With a probability Pi = 1/(L — 1) a tile from the gas 
hits site i, {i = 1, ...,L — 1). Depending on the value of the slope Sj at the site i, the 
following processes can occur: 

i) Si = and hi > hi_i. 

The tile hits a local peak and is reflected. 

ii) Si = and hi < hi_i. 

The tile hits a local minimum. With a probability Ma the tile is adsorbed 
{hi /ij + 2) and with a probability 1 — Ma the tile is reflected. 

iii) Si = l. 

With probability the tile is reflected after triggering the desorption of a layer of 
tiles from the segment {hj > hi = hi+b, j = i + 1, . . . ,i + b — 1), i.e. hj ^ hj — 2 
for j = i + 1, + 6—1. This layer contains 6—1 tiles (this is always an odd 
number). With a probability 1 — Wd, the tile is reflected and no desorption takes 
place. For an example see Fig. [91 

iv) Si = -1. 

With probability the tile is reflected after triggering the desorption of a layer 
of tiles belonging to the segment {hj > hi = hi^b, j = i — b + 1, . . . ,i — 1), i.e. 
hj v-^ hj — 2 for j = i — b + 1, ...,i — 1. With a probability l — u^ the tile is reflected 
and no desorption takes place. 

The physics of the model depends on one parameter u = 1/w = Ua/ud. We use 
these notations since we consider below either very small values of w (or u) and our 
intuition can be helped by the fact that for m = or w = the properties of the 
system are known. Namely, for both -u = and w = the spectrum is massive [1], the 
stationary states being the substrate {u = 0) and a full triangle {w = 0), i. e., hi = hi-i 
(z = l,...,L/2). 
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Figure 10. Phase diagram of tlie raise and peel model. For < u < 1 (phase I), the 
model is massive. At u Uc = 1 it is massless and conformal invariant. For 1 > w > 
it is scale invariant with varying critical exponents and exhibits SOC. For w ~ the 
system is massive. 



At the special point of the RPM where the ratio u = Ua/ud = 1/w = 1 (the 
Razumov-Stroganov point considered in earher sections) the model is governed by the 
Temperley-Lieb Hamiltonian (l28l) acting on the space of configurations by the ideal Ji 
defined in ( IT3il (see Sec. 4 and Appendix). As a consequence, the model is exactly 
integrable and conformally invariant. As discussed in Sec. 3 this last invariance allows 
exact predictions for the finite-size eigenspectrum of H. For u ^ 1 the model is not 
exact integrable anymore and our knowledge comes from from Monte Carlo simulations 



The phase diagram of the model obtained from the Monte Carlo simulations is given 
in Fig. [TOl In the domain < u < 1 the system is massive undergoing a second order 
phase transition at m = Mc = 1. For 1 > w > the conformal symmetry seen at w = 1 is 
broken but the system stays scale invariant. Since in this domain one has avalanches we 
call this phase a self-organized criticality (SOC) phase. Another phase transition occurs 
aX w = where the system becomes massive again. The phase diagram of Fig. [10] was 
established from the static and dynamic properties of the RPM [2]. In the next section 
we are going to consider the static properties of the stationary state and in Section 7 
the dynamical ones. 

6. The stationary state of the RPM 

It is useful to define some quantities which characterize the surface. An obvious 
geometric observable is the set of sites j {j even) for which hj = (the sites and 
L always belong to this set). These sites are also called contact points. This set is 
important, for example, to study desorption. Desorption events are limited to the area 
between two contact points, defined as a cluster (see Fig. [8]). The density of contact 
points gives a local observable for which one can define various correlation functions. 

We first define quantities in the stationary state. The average number of clusters 
k{L), the density of clusters p(L) and the local density of clusters g{j,L) are defined 



[2]. 



by: 



L L 



k{L) = Y.g{3,L) = {Y^5^J 



(37) 



p(L) = k{L)/L. 
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Figure 11. Typical configurations in the stationary states for L 
values of u. 



128 sites and three 



Note that the maximum value of p{L) is ^ which is obtained when u 



and one has 



only one configuration in the system: the substrate. The average height is defined as: 
h(L) = iS^. (39) 



The lowest value of h{L), corresponding to the substrate, is |. It is useful to consider 
also the average height in the middle of the system: 

h.iL) = {h.). (40) 

A relevant quantity is the average of the fraction of the interface where desorption does 
not take place (FND) 

L-l 



n{L) 



1 



L 



t(E(i-I^.I))- 

J- — 1 



(41) 



This quantity allows to estimate the average number of tiles desorbed in avalanches (see 
Sec. 7). By studying the behaviour of n{L) for large L and various values of w we will 
be able to establish the existence of a discontinuous phase transition at w = 0. 

The large L behavior of any of the quantities ( 137|) -(HT1). say h{L), will be denoted 



= lim h{L). 



(42) 



In the next section we are going to discuss average quantities which are not only 
L dependent but also time dependent. For these quantities we use, for example, the 
notation k{t, L) to denote the number of clusters. 

In order to get an intuition on the behavior of the model when u (or w = 1/u) 
changes, in Fig. [11] we show, for L = 128, typical configurations in the stationary states 
for three values of the parameter u. One notices that for u = 0.4 there are many clusters. 
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Figure 12. The density of clusters poo for various values of u < 1. 
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Figure 13. (a) The local density of contact points g{x, L) at m = 0.4 for different 
lattice sizes L {g{x,L) is zero for x even), (b) (Color online) The scaling function 
C{cos{-Kx/ L)) defined in (|46)) for u — 1. 



fewer for u = 1 and a single one for w = 0.4 (other, less probable, configurations do 
exhibit several clusters). We also notice that the average heights are quite small. 

We first consider the density of clusters p{L) in the domain < m < 1. Taking 
large values of L, one obtains the u dependence of poo shown in Fig. [121 In all this 
domain the density of clusters stays finite. The density of clusters decreases from its 
maximum possible value 1/2 which corresponds to the substrate at u = to the value 
zero for u = 1. This indicates that there is a phase transition at m = 1. From the Monte 
Carlo simulations it was hard to obtain the critical exponent which gives p{u) when u 
approaches the value one. 

The local density of contact points g{x,L) (see Eq. ( l37j) ) for u = 0.4 and different 
lattice sizes is shown in Fig. [TST a). One observes that for small values of x, and large 
values of L, g{x,L) decreases exponentially from the value 1 at x = (not shown on 
Fig. fT3l(a)) to a constant value in the bulk. This value is twice that seen in Fig. [121 
since for even x one always have g{x) = 0. 

Since for u = one can diagonalize the related Hamiltonian of the stochastic process 
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and show rigorously [IJ that the system is massive, we can expect that the system is 
massive in the whole domain < u < 1. 

We turn now to the interesting case u = 1 where the system is exact integrable and 
conformally invariant. In Ref. [8j a conjecture was made which gives the probability 
Pn{k) to have k clusters in a system of size L = 2n: 

^ 4-+^ r(i + n + A:)r(i)r(3n + 2)r(2n-A:) 

ny ) ^ 27- r(i)r(| + 2n + k)V{n + l)V{2n + k + 2)V{n - k + 1)'^ ' 
The last result gives us the leading behaviour 

k{L) = ±kP^{k)^^-^^I^Ll, (44) 
k=i 27r 

for the average number of clusters. This shows that the density of clusters vanishes 
algebraically for u = 1. This conjecture was checked using Monte Carlo simulations 
up to a system of size L = 512. Since, as we show below, the density of clusters also 
vanishes algebraically in the entire domain i;; < 1, it is convenient to define a critical 
exponent a < 1 which gives the power law increase of the number of clusters with the 
system size, k{L) ~ L". For -u = 1 we have obviously a = 2/3. 

It is interesting to consider, again for u = 1, the density of contact points g{x, L) 
for several lattice sizes and notice that for the function 

G{x/L) = L\g{x,L), (45) 

one has a data collapse. As suggested by the conformal invariance for the confined 
critical systems [33] we make an ansafe about the functional behaviour of G{x/ L) taking: 

g{x, L) = [- sm{TTx/L)]-^/^C{cos{7ix/L)). (46) 

In Fig. [T^ b) one shows the function C(cos(7rx/L). It is practically a constant. If one 
assume that C is independent of ttx/L, one can get its value by integrating g{x,L) 
(Eq. (H6ll ) and compare the result with the average number of clusters given by (^^. 
One obtains: 

C = -^r(-^) = 0.753149..., (47) 
ovre t) 

in very good agreement with the data obtained from Monte Carlo simulations (see 
Fig.IHb)). 

The motivations for the ansatz fj46l) comes from the following observation. If the 
density of contact points corresponds to a local operator in conformal field theory, one 
expects [33] . 

g{x, L) = C[- sm{7ix/L)f, (48) 
vr 

where X = A + A is the scaling dimension of the local operator (A — A is the 
conformal spin). Comparing Eqs. fHBl) - fH51) one concludes that the density of contact 
points corresponds to a local operator with scaling dimensions X = 1/3. This result 
is surprising. One can show that the profile of an operator with spin vanishes 
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w 


1 


0.85 


0.75 


0.4 


0.25 


0.1 


0.05 


0.025 


a 


0.67 


0.50 


0.40 


0.30 


0.24 


0.15 


0.06 


0.01 



Table 1. Estimates of the exponent a giving the increase of the number of clusters 
with the size of the system L. These data were obtained from numerical analysed of 
several lattice sizes up to L = 65536. 

This implies A = A = 1/6. This value can't be obtained from the possible values of 
A (see fl20l) ) of minimal models. It might be obtained from algebras [Ml [35] 

but it is not clear why W-symmetry should play any role in our problem. Another 
possible explanation is that for a conformal theory with c = (this is also the case 
of percolation), under certain circumstances, the scaling dimension can be an arbitrary 
number [361 [37]. This possibility is a bit strange since the number one has to obtain is a 
neat 1/6. Finally, it is possible, that the proof that operators with spin have no profile 
functions does not apply to our problem. If this is the case, one can choose A = and 
A = 1/3 The value 1/3 is a perfectly acceptable one (see Eq. (123|) in Sec. 3). To 
conclude, we have not yet an explanation for the value X = 1/3 in Eq. (HHl) . 

We move away from u = 1 and consider the domain 1 > w. The number of clusters 
k{L) keep having an algebraic increase with an exponent a that decreases monotonically 
with w. In Table [1] we give estimates for the exponent a for several values of w. It was 
also observed that for several values of w the scaling law 

G{x/L) = L^~''g{x,L) (49) 

holds. 

These results indicate that for m > 1 the system stays scale invariant with a varying 
critical exponent. 

Although it seems obvious that one should look at the L dependence of the average 
value of the heights h{L) or the average value of the height at the middle of the lattice 
hi{L), it turns out that these quantities are harder to study since they vary slowly with 
the system size (see [2] for some illustrating figures). In the case of hi{L) we have the 
dependence shown in Fig. [TU^a). 

A different picture emerges if one considers very small values of w as shown in 
Fig. [T^ b). One can observe two different phenomena. The values of hi{L) increase 
substantially as compared to the values observed at larger values of w, and more 
interestingly, it looks like they saturate for large values of L. We believe [2] that this 
saturation phenomena happens in the whole region u > 1 {w < 1). We did not see such 
saturation for smaller values of u (see Fig. [T4[(a)) because the lattice sizes considered in 
the simulations was not large enough for these values of u. Considering these simulations 
with lattice size L = 65536 a good fit, for various values of if, is obtained by the curve 

§ For the appearance of operators with conformal spin in the presence of the quantum group symmetry 
mentioned in Sec. 3, see Ref. [38] , 
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Figure 14. (a) Average height in the middle of the lattice hi{L) as a function of 
the lattice size L for some values of w in the range 1 > w > 0.025. (b) Average height 
in the middle of the lattice hi{L) as a function of the lattice size L for various small 
values of w (w < 1/200). 



m 



4.67 



w 



0.62 



(50) 



Interestingly, if we use ( l50l) for w = 1 we obtain /imax = 4.67, a value compatible with 
the results shown in Fig. [MT a). Taking the results shown in Fig. [T4l(b) at face value, 
implies that the phase SOC (see Fig. [101) extends at least to the value w = 10~^. On 
the other hand for w = the single configuration is the full triangle, /imax has values 
of the size of the system and therefore the L-dependence of h^ax shows no saturation. 
This suggests a non-analytical behaviour of various quantities describing the system at 
w = 0. In order to illustrate this point, in Fig. [15] we show typical configurations in 
the stationary state for a very small value of w [w = 0.00025) and different system 
sizes. One can see a clear change in the role of the boundaries if L varies from 1024 to 
32768. For L = 1024 one is close to the full triangle configuration (expected for w = 0). 
However, as one increases L, one gets a plateau (in which the height is almost constant 
and L independent) and the boundaries play a less and less important role. 



7. Space-time phenomena and avalanches 

In the last section we present results that corroborate the phase diagram shown in 
Fig. [TOl In the present section we give additional facts supporting these results. We are 
going to use the Family- Vicsek [39] scaling in order to check if the system is in a scale 
invariant phase and if the answer is positive, determine the dynamical critical exponent 
z. If we consider a time-dependent average quantity a{t,L), being a(L) the average in 
the stationary state, than for large values of t and L the function A{t, L) scales as: 

a(t,L) . t 



A{t,L) 



1 ~ A{- 



(51) 



a{L) ' "'L' 

where z is the dynamical critical exponent. The scaling function depends on the initial 
conditions. In the data shown below, we have taken as initial condition the substrate 
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Figure 16. (Color online) (a) Function K{t,L), defined on (|5ip . for u ~ 0.4 as a 
function of time for the lattice sizes L = 2048, 4096 and 8192. (b) Function K{t, L) for 
M = 1 as a function oit/L for various lattice sizes. Red curve is given by a fit to ([?T|) 
considering only the first exponential with Ai — ttSVS and Ci — 0.18. In the inset we 
compare K{t, I) - Cie"^^*/^ for L = 8192 with the short time exponential C2e~^^*/^, 
given in red line (color online) with A2 = |Ai and C2 = 0.54. 



[h{2i) = 0,h{2i + 1) = 1, {i = 0, ...,n - 1), Hl = 0) with probability one. We 
first consider the quantity K{t, L) which corresponds to the average number of clusters 
{a{t,L) in (ISTl) is k{t,L) in this case). 

In Fig. [T6](a) we show the function K{t, L) for u = 0.4 for different lattice sizes. We 
notice that for large lattice sizes, K{t, L) is a function on t only and can be obviously 
fitted with a sum of exponential functions as expected in a massive phase. Fig. [T6](b) 
shows K{t,L) for u = 1 and as expected from conformal invariance, z = 1. Conformal 
invariance gives also information on the function K{t, L). One can use the knowledge of 
the finite-size scaling limit spectrum and obtain the functional dependence of K{t/L). 
To illustrate how it works see Fig. [T6l(b). For large values of t/L one can do a fit to the 
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Q7 


0.9 


0.86 


0.74 


0.72 


0.7 


0.55 


0.52 


0.50 


0.4 


0.35 


0.34 


0.34 


0.25 


0.31 


0.30 


0.30 


0.025 


* 


0.07 


0.07 



Table 2. Estimates of the dynamical critical exponent z for various values of w. 
The estimates were obtained using ([51]) with A replaced by K, H and Hi. The 
number represented by (*) is not reliable since in that region the number of clusters 
is very small. These estimates were obtained by considering several lattice sizes up to 
L = 65536. 



Alt Aot 



data, using the following ansatz: 

K{t/L) = Cie-^ + Cae-^, (52) 
where the Aj are obtained from the spectrum of the Hamiltonian (l28ll . given by f l22l) : 

E, = ^, 1 = 1,2. (53) 

The constants Ci depend on the initial conditions. If one starts with the substrate 
Ci = 0.28 and C2 = 0.54. 

We now consider smaller values of w (larger values of u). We have done an extensive 
study of various functions A(t,L) taking K{t,L), H(t,L) and Hi(t,L) in (ISTj) . The 
estimated values for the critical exponent z are shown in Table El For very small values 
of w, the estimates coming from K{t, L) are poor because the number of clusters is 
small and hence large errors. Inspection of table [2] shows a smooth drop of the values 
of z from z = 1 for w = 1 to 2; ~ for w ^ 0. The latter values is in agreement 
with a direct calculation of the mass gap at w = [T] which gives a finite gap in the 
thermodynamic limit. In the region where z —>■ the information propagates, through 
the system, at infinite velocity. This can be explained nicely by Fig. [15] for large L. If 
the dynamics chooses a site near the edges, the peel process affects almost the whole 
interface, simultaneously lowering the entire upper plateau by a unit, implying that the 
information travelled an infinitely long distance in the limit L ^ 00 |40j . 

We have not studied extensively the two-contact points space-time correlation 
functions c{R,t,L), except for the special case u = 1. This case is interesting since 
from conformal invariance it exhibits, far from the boundaries, the following scaling 
form in the bulk (see ([15])): 



2 



C{R,t,L) = L--^G{iJ.), (54) 
where 



^^ = ^ (55) 
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Figure 17. (Color online) The two-contact point correlation function C{R, t, L) times 
Lt as a function of the scaled distance R/ L and the scaled time {vg — |\/3 is the 
sound velocity) for various lattice sizes. 



In f l5^ and fl55l) R is the distance between the contact points, t is the time difference 
and Vs = |-\/3 is the sound velocity [2]. The factor L~3 in (IMl) is obtained from ( H5l) . 
In Fig. [T7| we show the scaling function G{fi), where we have used the unconnected 
correlation function in (JSj). To avoid boundary effects the data was taken only in the 
center segment of size L/2. 

The functional dependence (|5^ which is special for the conformal invariant case is 
nicely seen in Fig. [T71 Our data are not good enough in order to extract from the small 
H behavior of G{ii), the two-contact point correlation function in the thermodynamic 
limit. 

An interesting property of the RPM is the occurrence of avalanches in the desorption 
processes. For the case u = 1, based on exact results on small lattices, the existence 
of avalanches was suggested in [1]. The Monte Carlo simulations on large lattices show 
avalanches in the whole domain < w; < 1. 

In the stationary state, once a tile from the rarefied gas hits the interface, it can be 
reflected, adsorbed or can trigger a nonlocal desorption process in which many tiles may 
leave the interface, this defines an avalanche. In the latter case the size of the avalanche 
is given by the number of tiles T that are released in the process. This number is always 
odd, therefore it is convenient to write: T = 2v ~ 1 {v = 1,2, . . .). For m < 1 only a 
finite number of tiles are removed, since the density of clusters is finite. For w < 1, the 
cluster density vanishes, and therefore macroscopic number of tiles may be desorbed, 
hence macroscopic avalanches. 

We denote by S{v, L) the PDF which gives the probability for an avalanche of size 
V for a system of size L. In models of self-organized criticality (SOC) [4Tj-[45j, for large 
values of v and L, one may expect the PDF to exhibit a simple scaling law, 

S{v,L)=v-^F{^), (56) 
characterized by two exponents r and D. 
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Figure 18. (Color online) The scaling function F{v/L) for w = 1. The data are for 
lattice sizes L = 1024 and L = 2048. 
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0.95 


0.90 
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0.30 


0.05 


0.02 


0.005 


D 


1.004 


1.089 


1.087 


1.066 


1.040 


1.026 


1.017 


1.002 


1.002 


1.006 


T 


3.00 


2.77 


2.63 


2.47 


2.32 


2.25 


2.18 


2.07 


2.046 


2.00 



Table 3. Estimates of the critical exponents D and r obtained by using lattices sizes 
L = 4096 and L' = 8196. 



The Monte Carlo simulations |2j allow the evaluation of the above exponents. For 
M = 1 we have r ~ 3 and D 1. The result -D = 1 is to be expected since in a conformal 
invariant theory (this is the case if u = 1) one has no other length than the size of the 
system L. We found no explanation for the value three of r although this number shows 
up in many aspects of the model (see for example (l45l) ). In Fig. [18] we show the scaling 
function F{v/L) defined in fl56l) . One observe a nice data collapse. 

In Table |3] we give the estimates for the exponent r and D for various values of w. 
The estimates presented in this Table are obtained using only one pair of sizes (4096 
and 8192) but we made sure that the results are reliable by studying the changes of the 
values of the estimates using smaller lattices. 

The results of Table [3] indicate that the exponent D stays unchanged and equals 
to one in the whole range w <1. This can be understood in the following way [16] the 
tiles which are desorbed and create the avalanche belong always to a one- dimensional 
layer. 

We turn now to the exponent r. As can be seen in Table [3] this exponent varies 
from the value 3 (m = 1) to a value close to 2 for w close to zero. In the case u = 1, the 
average size of the avalanches stays finite when L gets large but the dispersion diverges 
(the function F{v/L) shown in Fig. [18] does not vanish at the origin). Does an exponent 
r = 2 mean that the average size of the avalanches diverges logarithmically with L ? Not 
necessarily, it depends if the scaling function F[v/L) vanishes or not at the origin. In 
Fig. [I9]we show for w = 0.00025, that F{v/L) does not vanish at the origin. Since later 
in this section we are going to show that {v) stays finite for any w > 0, we conclude that 
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Figure 19. (Color online) The functions F{v/L), defined on ((56| . a.t w — 0.00025 for 
some values of L. 



T gets very close to the value 2 but never reaches it. For a more precise determination 
of the exponent r, see Figs. 15-18 in Ref. [2]. 

The Monte Carlo simulations show avalanches at least in the interval 1 > w > 
0.00025 and therefore the system is scale invariant in this interval. It is probably safe 
to assume that the SOC phase covers the domain 1 > w > 0. It was checked that for 
u < 1, where the system is massive, for large values of v the PDF has an exponential 
decrease and not an algebraic one. 

We turn now our attention to the average size of the avalanches. A simple mean- 
field argument allows to compute it using informations about the stationary state only. 

We first consider, in the stationary state, the average fraction of interface where 
desorption does not take place n{L) (FND) defined in (HTl) . This quantity varies from 
the value 1 for the substrate (the single configuration for m = 0) to the value 1/(L — 1) 
for the full triangle (the single configuration for w = 0). 

For u = 1 there is a conjecture for the values of n{L) 



This conjecture was checked using our Monte Carlo simulations for various lattice sizes. 

If one knows n(L), the average probability Pa{L) to have an adsorption process, 
the average probability Pr{L) to have a reflection process and the average probability 
Pd{L) to have a desorption process can be easily obtained: 




2L + 2 



3 



1 



(57) 




n(L) 



1 




2 





(59) 



P,{L) + Pr{L) + P,{L) = 1, 



(60) 



where we choose Ua = I and Ud = w = 1/u. 
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Taking into account that in the stationary state the average number of adsorbed 
tiles is equal to the average number of desorbed tiles, one obtains: 



Note that this is a mean-field calculation since we have first computed the average 
probability to have a desorption process and multiplied it with the average number of 
tiles which are desorbed. In the large L limit (!58!) - fl6T!) gives: 



In Fig. [20] one shows the w dependence of as obtained from the extrapolated 
results of our Monte Carlo simulations. One notices a discontinuous behavior of n^o 
around w = 0. At w = where in the stationary state one has only one configuration, 
which is the full triangle (see Fig. IA2I) . one has = 0, while lim^^,^o+ ^oo = 0.5. Using 
f l62p one concludes, in agreement with previous observations, that for any finite value 
w, the average size of the avalanches is finite. We consider that the discontinuity of 
at = shown in Fig. [20] is one of the most interesting and unexpected property 
of the model. The value noo = 0.5 observed for w = 0+ corresponds to an interface in 
which on the average we have an equal number of sites where one has local minima, 
maxima, positive or negative slopes. It is not clear to us why this interface is special. 
The existence of this discontinuity also gives supplementary support for the existence of 
the SOC phase. Studying the finite-size scaling properties of the Hamiltonian spectra 
for small lattices and different values of w in Ref. [1] it was shown that for w < 1 one 
does not have conformal invariance. On the other hand, for w = 0, we know [Ij that the 
system has an energy gap. In spite of all our determinations of critical exponents in the 
domain < w < 1 implying that the system is gapless, one could argue that one sees 
crossover effects and that one can have a massive phase in the whole domain < w < 1. 
The discontinuity observed in the behaviour of noo assures us that this is not the case. 

In Table [1] we compare in two cases the mean-field predictions for (T) l given in 
f [62|) with the values measured in Monte-Carlo simulations. The agreement is excellent. 
In the same Table the values for n(L) for several values of L are also given. One notices 
that for M = 1 and L = 1024 one already gets, within five digits, the asymptotic value 
3/4. This is not the case for w = 0.00025, where the finite-size effects are larger. Our 
numerical analysis indicate that the asymptotic behavior of n(L) can be described by 



where ^4(1) = 1/6 and A{w) diverges if w —>■ (data not shown). 

We conclude this section with the following observation, which is necessary for the 
consistency of the present model for avalanches. If one looks at the expression ([59]) of 
the average (over configurations) probability that a tile hitting the interface is reflected, 
one sees that for small values of w, the probability approaches the value one. This 
implies that to obtain large avalanches (small r) one has to wait a long time. This is to 
be expected in models of SOC. 




(61) 




(63) 
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Figure 20. The w dependence of the fraction of the interface where adsorption does 
not take place Uoo obtained from the extrapolated results of Monte Carlo simulations. 



w 


L 




{T)l 


MF 


{T)l/MF 


1 


1024 


0.749 


1.497 


1.497 


0.999 


1 


2048 


0.749 


1.498 


1.501 


0.998 


1 


4096 


0.749 


1.499 


1.500 


0.998 


1 


8192 


0.749 


1.499 


1.499 


0.999 


0.00025 


1024 


0.199 


496.999 


496.119 


1.001 


0.00025 


2048 


0.313 


914.346 


913.820 


1.001 


0.00025 


4096 


0.407 


1374.713 


1373.856 


1.001 


0.00025 


8192 


0.455 


1674.877 


1674.550 


1.000 



Table 4. Values of and the average number of tiles obtained from Monte Carlo 
simulations {T)l as compared with the mean- field results (MF). The values are given 
for several lattice sizes L and two values of w: w = 1 and w = 0.00025. 



8. Conclusions 

We have shown that the one parameter u = 1/w dependent RPM has an unusual phase 
diagram (see Fig. [TOj) . At -u = 1 (the Razumov-Stroganov point) the model has many 
facets. Firstly, the stationary state is related to a two-dimensional equilibrium system. 
The details of this observation are still at the level of a conjecture. The fact that 
the normalization factor of the stationary PDF is related to an enumeration problem 
of a two-dimensional system (fully packed loops or alternating sign matrices) was 
recently proved rigorously [IE]- [50]. The two-dimensional system has many interesting 
combinatorial properties which are reflected in the PDF of the stationary state. Also 
at u = 1 the Markov process has a space-time symmetry (conformal invariance) which 
does not exist for u 1. 

How does the combinatorial facet of the model can be used to derive meaningfuU 
physical quantities? Using small system sizes, one can "guess" expressions which can 
be checked by Monte-Carlo for large systems. One example is the probability to have k 
clusters for a system of size L = 2n (see Eq. (HI])). Actually this expression is a solution 
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of some bilinear recurrence relations (Pascal hexagon [IB])- For another example, see 
Eq. dSH]). 

The conformal invariance facet of the model has other consequence. It fixes the 
value of the dynamical critical exponent: z = 1. It also fixes the finite-size scaling 
behavior of the density of contact points (see Eqs. (H71) and (H8l) ). the form of the two- 
contact point correlation function (see Fig. [T71) and the Family- Vicsek scaling functions 
(see Eq. ([53]) and Fig. [Te^b) ). 

The RPM has been extended in several directions. In one of them, one has kept 
the Hamiltonian (!28|) but enlarged the vector space on which it acts to a larger left 
ideal. This has allowed us [51], keeping conformal invariance, to study the move of a 
random walker in an unquenched disordered system as well as the annihilation reaction 
A + A — s> in the same medium. The same study for other models, without conformal 
invariance, can be found in [52] and [53] . 

Some very interesting results have been obtained by extending the TL semigroup 
algebra to the one-boundary TL semigroup algebra [I7]-[19]. The configuration space is 
different and biased toward the one end of the system where a new generator is placed. 
The model gets new combinatorial facets and conformal invariance is maintained. 

Before closing this paper, we would like to comment on the phase diagram of the 
RPM (see Fig. [TOj) . If we perturb the system and go away from u = 1, one obtains 
a massive phase if u < 1. This happens with many conformal invariant models. If 
however we take u > 1, the system stays scale invariant with varying critical exponents 
without being conformal invariant. This phenomenon is new and it is a consequence of 
the nonlocality of the model [51] . From a phenomenological point of view we think that 
the observation that for any finite u the average height stays finite for large lattice sizes 
is very unexpected (for unweighted RSOS paths it increase like L^/^). Keeping in mind 
that for w = 0, in the stationary state, there exists only one configuration which is a 
full triangle (see Fig. IA2I) one is tempted to believe that configurations where on "feels" 
the triangle effects should play an important role. This is not the case (see Fig. [T5|l . 
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Appendix. Link diagrams and RSOS (Dick) paths 

We consider a one-dimensional lattice with L + 1 [L = 2n) sites. The configuration space 
of the fluctuating interface is given by RSOS paths (better known in combinatorics as 
Dick paths). They are defined by taking nonnegative heights hi which obey the restricted 
solid-on-solid (RSOS) rules 

K+i -hi = ±1, ho = hL = 0, hi > 0. (A.l) 

See Fig. IA2I for an example. 
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Each RSOS path corresponds to a hnk pattern (see section 3) in the following 
way. On each RSOS path, draw the equal height contour lines as in Fig. lAlb . By 
straightening out the surface, keeping the contour lines and rotating the picture around 
the horizontal axis, we obtain Fig. lAlb . The contour lines connect pairs of sites and 
Fig. lAlb defines a link pattern. Notice than for an RSOS path defined on L + 1 sites 
corresponds a link pattern defined om L sites. 

(a) ^^^^^^^..^ry^^r:^=^r:;F^ (6) 

Figure Al. An interface with contour lines (a) and the corresponding link pattern 
(b). 

In the case L = 6 for example, to the five link patterns (129|) correspond five RSOS 
configurations 





Figure A2. The RSOS configurations corresponding to the links patterns of the 
L = 6 lattice {{L + 1) points). 



One can read off the action of the Hamiltonian (1301) on this configuration space and 
therefore get the different transition rates. These are the transition rates of the RPM 
at the Razumov-Stroganov point used in Section 5. 
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